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§ 1 .  Introduction . 


This  report  presents  an  analysis  of  the  structure  of  transient 
acoustic  waves  in  stratified  fluids  whose  densities  and  sound  speeds  are 
functions  of  a  single  depth  coordinate.  The  acoustic  field  in  such  a 
fluid  may  be  described  by  a  real-valued  function  u(t,x,y)  (the  acoustic 
potential  or  the  excess  pressure)  that  satisfies  the  wave  equation 

(1.1)  It^  ~  c^(y)  p(y)  ^  *  (p~'(y)7u)  -  f(t,x,y) 

where  t  is  a  time  coordinate,  y  is  a  depth  coordinate,  x  (X),X2)  are 
Cartesian  coordinates  in  a  horizontal  plane,  7  *  (3/3xi ,3/3x2, 3/3y)  and 
f(t,x,y)  is  a  function  that  characterizes  the  wave  sources.  c(y)  and 
p(y)  are  the  variable  sound  speed  and  density,  respectively,  and 
P~^ (y)  •  l/P(y) . 

This  report  is  a  sequel  to  the  author's  report  [12]  on  "Spectral 
Analysis  of  Sound  Propagation  in  Stratified  Fluids."  That  work  contains 
a  spectral  analysis  of  the  acoustic  propagator 

(1.2)  Au  -  -c^(y)  p(y)  7  ♦  (p  \y)  7  u) 

for  the  cases  of  an  unlimited  fluid  (-«  <  Xj.x^.y  <<*>),  a  semi-infinite 
layer  (-“  <  <  «  and  0  <  y  <  »)  and  a  finite  layer  (-<»  <  XjjX^  <  ® 

and  0  <  y  <  h  <  <») .  The  Integration  of  (1.1)  below  is  based  on  the 
spectral -analysis  of  [12]  and  the  notation  and  results  of  [12]  are  used 
throughout  this  report.  As  in  [12],  only  the  case  of  an  unlimited  fluid 
is  presented  in  detail.  The  modifications  required  in  the  second  and 
third  cases  are  described  in  §6  at  the  end  of  the  report. 
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The  analysis  of  the  structure  of  transient  acoustic  waves  in 
plane  stratified  media  was  initiated  by  the  author  in  1973  [7].  This 
work  dealt  with  the  special  case  of  the  Pekeris  profile  (propagation  in 
a  half-space,  p(y)  ■  const.,  c(y)  piecewise  constant).  In  1974  the 
results  were  extended  to  the  symmetric  Epstein  profile  (p(y)  ■  const., 
c  ^(y)  “  K  sech*  (y/H)  +  M)  [8].  The  general  Epstein  profile 
(p(y)  “  const.,  c  ^(y)  “  K  sech^  (y/H)  +  L  tanh  (y/H)  +  M)  was  treated 
in  1979  [11]  using  the  spectral  analysis  of  the  Epstein  operator  due  to 
Guillot  and  Wilcox  (2,  3].  A  preliminary  version  of  [11]  was  announced 
in  1978  [I].  The  analysis  of  these  cases  was  based  on  explicit 
representations  of  their  normal  mode  functions  by  means  of  well-known 
special  functions.  Extension  of  the  analysis  to  larger  classes  of  strat¬ 
ified  fluids  had  to  await  the  extension  of  the  spectral  analysis  of  the 
acoustic  propagators  (1.2)  to  such  classes. 

The  class  of  stratified  fluids  treated  in  [12]  was  characterized 
by  the  properties 


(1.3) 


(1.4) 


(1.5) 


p(y)  and  c(y)  are  Lebesgue  measurable, 

0  <  <  P(y)  <  Pm  <  «.  0  <  c^  <  c(y)  <  Cm  <  », 


p(±<»)  I 


dy  <  <*>, 


•+0O 

|c(y)  -  c(±«>)  I  dy  <  », 
JO 


where  p^,  p^,  p(±*),  Cj^^,  Cm  and  c(±<»)  are  constants.  In  this  report  the 
results  of  [10]  and  [11]  on  transient  acoustic  waves  in  stratified  fluids 
are  extended  to  the  class  of  fluids  characterized  by  (1.3),  (1.4)  and 
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(1.6) 


lp(y)  -  P(±«>)  1  <  c(±y) 


. 

lc(y)  -  c(±<»)l  <  C(±y)"“ 


for  ±y  >  0, 


where  C  and  a  are  constants  and 


(1.7)  a  >  I  . 

It  Is  clear  that  (1.6)  Implies  (1.5).  It  will  be  seen  from  the  analysis 
below  that  (1.6)  could  be  replaced  by  other  order  conditions  at  y  = 

It  is  not  known  whether  the  results  of  this  report  hold  for  the  entire 
class  of  lulds  defined  by  (1.3),  (1.4),  (1.5). 

The  remainder  of  this  report  is  organized  as  follows.  In  §2  the 
normal  mode  expansions  for  unlimited  fluids  of  [12]  are  used  to  decompose 
transient  acoustic  fields  into  free  and  guided  components  and  to  obtain 
integral  representations  of  these  components.  The  behavior  for  large  t 
of  the  free  and  guided  components  is  derived  in  §3  and  §4,  respectively. 
In  §5  these  results  are  used  to  calculate  the  asjmptotlc  distribution  for 
large  t  of  the  wave  energy.  §6  presents  corresponding  results  for 
semi-infinite  and  finite  layers. 


§2.  Normal  Mode  Expansions  of  Transient  Acoustic  Fields. 

The  integration  of  the  wave  equation  (1.1)  will  be  based  on  the 
spectral  analysis  of  the  acoustic  propagator  A  in  the  Hilbert  space 
3C  ■  LjCR^iC  ^(y)  p  ^(y)  dxdy) ,  as  developed  in  [12].  To  this  end  (1.1) 
is  interpreted  as  the  equation 

(2.1)  +  Au  =  f(t,*),  t  6  R, 


for  a  function  t  -*■  u(t,*)  e  3C.  The  wave  sources  will  be  assumed  to  act 
during  a  time  Interval  [0,T],  so  that  supp  f  C  [0,T].  The  corresponding 
acoustic  wave  is  the  solution  of  (2.1)  that  satisfies  the  initial 
condition 

(2.2)  u(t,*)  »  0  for  all  t  <  0. 

The  solution,  based  on  the  spectral  theorem  for  the  selfadjoint  operator 
A  >  0,  is  given  by  Duhamel's  integral 

(2.3)  u(t,*)  -  I  {A  sin  (t  -  t)a‘''^}  f(T,*)dT,  t  >  0. 

Indeed,  if  f  €  C([0,T],3f)  then  (2.3)  is  the  unique  "solution  with  finite 
energy"  of  [6],  while  if  f  6  C([0,T],D(A^/^  ))  then  (2.3)  is  the  "strict 
solution  with  finite  energy."  In  addition,  if  f  6  C( [0,T] ,D(A~^/^  ))  then 

(2.4)  u(t,x,y)  »  Re  {v(t,x,y)} 
where  v(t,*)  is  the  complex-valued  potential  defined  by 
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(2.5)  v(t,‘)  »  i  exp  |  exp  {ItA^^^}  f(T,*)dT. 


In  particular. 


(2.6) 


v(t,*)  ■  exp  {-itA*^^}  h  for  all  t  >  T 


where 


(2.7) 


«p 

i  A  exp  {ItA^^^}  f(T,*)  dt. 


The  Initial  value  problem 


(2.8) 


+  Au 


0  for  t  >  0, 


(2.9) 


u(0) 


f  du(0) 
’  dt 


g 


can  be  treated  by  the  same  formalism.  Indeed,  if  f  £  D(A^^^  )  and 
g  €  D(a”^^*  )  then  the  solution  of  (2.8),  (2.9)  is  given  by  (2.4),  (2.6) 
with  h  -  f  +  i  a"^^  g  €  D(A‘/2  j.  ^9^  d,.  3]. 

The  integral 

(2.10)  E(u,K,t)  -  I  {ivu|^  [1^]  ^ 

may  be  Interpreted  as  the  energy  of  the  acoustic  field  u  in  the  set 
K  C  R*  at  time  t.  Moreover,  A  is  the  selfadjoint  operator  in  X  asso¬ 
ciated  with  the  sesquilinear  form  A  on  3C  defined  by  D(A)  ■  1^(8^)  C  Jf 


and 
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(2.11) 


A(u,v) 


Vu  •  7v  p  ^ (y)  dxdy 

4* 


It  follows  from  Kato's  second  representation  theorem  [4,  p.  331]  that 
)  =  L^(R^)  and  for  all  u  €  ) 


(2.12)  ul^  =  A(u,u)  =  I  lVu|^  p  ^  (y)  dxdy. 

Hence  the  total  energy  satisfies 


(2.13) 


E(u,R^t)  =  liA^2 


3u 

3t 


If  h  €  D(A^''^)  and  u  is  defined  by  (2.4),  (2.6)  then  a  simple  calculation 
shows  that 

(2.14)  E(u,R^,t)  *  OA^^^  hll^  for  all  t  >  T. 

The  analysis  of  the  structure  of  the  acoustic  potential  (2.4), 
(2.6)  presented  below  is  based  on  the  normal  mode  expansion  for  A  of 
[12,  §8].  The  orthogonal  projections  {P^,P_,P5,Pi ,P2, • • •}  in  3f  defined 
by  the  normal  modes  form  a  complete  family  that  reduces  A  [12, 
Corollaries  8.4  and  8.10].  Hence  the  same  is  true  of  the  orthogonal 
projections 


(2.15) 


Pf  “  P+  +  P_  +  Po 


and 


P 

g 


No-l 

I 

k*l 


(2.16) 


This  provides  a  decomposition 


(2.17)  u(t,*)  =  u^(t,')  +  Ug(t,*) 
into  orthogonal  partial  waves 

(2.18)  Uj(t,*)  =  Pj  u(t,*).  Ug(t,*)  =  Pg  u(t,*). 

The  function  u^,  called  the  free  component  of  u,  will  be  shown  to  behave 
for  large  times  like  a  diverging  spherical  wave.  The  function  u  ,  called 

O 

the  guided  component  of  u,  will  be  shown  to  be  approximately  localized 
near  one  or  more  of  the  planes  y  *  const,  where  c(y)  has  local  minima  and 
to  propagate  outward  in  horizontal  planes  like  a  diverging  cylindrical 
wave.  This  second  component  shows  the  profound  effect  of  acoustic  ducts 
on  transient  acoustic  waves.  It  is  absent  when  c(y)  has  no  local 


minima.. 


§3.  Transient  Free  Waves. 


The  normal  mode  expansions  of  [12]  are  used  in  this  section  to 
calculate  the  asymptotic  behavior  for  t  “  of  the  free  component 
u^(t,*)  *  ?£  u(t,*).  The  principal  result  is  that  in  each  of  the  half¬ 
spaces  R^(d)  and  R^Cd),  where  R|(d)  -  {(x,y)  :  t(y  -  d)  >  O},  Uj(t,*)  is 
as}nnptotically  equal  to  a  wave  function  for  a  homogeneous  fluid  with 
propagation  speed  c(“>)  and  c(-“>),  respectively.  It  is  this  behavior 
that  motivates  the  term  "free  component"  for  u^Ct,*). 

It  will  be  asstjffled  that  the  total  acoustic  potential  u  satisfies 

u(t,*)  •  Re  {v(t,*)}  where  v(t,*)  ■  exp  {-itA*^^}h  and  he  D(A^^^) 

(see  §2).  The  corresponding  partial  waves  Uj^(t,*)  ■  u(t,*)  with  k  >  1 

satisfy  Uj^(t,*)  -  Re  {exp  (-it  A*^^)  This  follows  from  the  fact 

that  the  normal  mode  functions  <j'jj(y»p)  are  real  for  k  >  1,  which  implies 

that  Pj^(h)  =  follows  by  addition  that  Ug(t,*)  * 

Re  {exp  (-itA^^^)  P  h}  »  Re  {v„(t,*)}  and  hence 
S  g 

[  Uj(t,*)  -  Re  {v^(t,*)}. 


(3.1) 


v^(t,*)  -  exp  (-it  A*''^)  P^h  -  exp  (-it  A^^^)  h^. 


The  normal  mode  representations 


(3.2)  v,(t,x,y)  -  1  (f+Cx.y.P.q)  exp  {-ItX'^*  (p,q)}  h+(p,q)  dpdq 

r  )r3  - 

of  [12,  §8]  provide  the  starting  point  for  calculating  the  asymptotic 
behavior  of  Uj(t,x, y)  for  large  t.  The  integral  in  (3.2)  will,  in 
general,  converge  only  in  X.  For  brevity  the  3C-lim  notation  of  [12]  is 
suppressed  in  this  report. 
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Equation  (3.2)  gives  two  representations  of  corresponding  to 
the  two  families  and  The  calculations  below  are  based  on  the 

(|>_-representation  which  has  been  found  to  yield  the  simplest  form  of 

the  asymptotic  wave  function.  It  will  be  convenient  to  Introduce  the 
characteristic  functions  x^,  Xo  )(_  cones  C^,  Cq  and  C_  in 

(p.<l)~space  [12,  (1.48)-(1.51)  ]  and  to  decompose  h_  as 

(3.3)  h_(p,q)  =  2.(p,q)  +  m(p,q)  +  n(p,q) 

where  ■  )Cj.  h_,  m  *  Xo  “  “  )(_  ^_*  corresponding  decomposi¬ 

tion  of  v£  is 

(3.4)  ^f  -  +  ^m  +  ^n 

where 

' 

Vj^  =  exp  (-it  A*^^)  $*Jl 

(3.5)  ■  V  *  exp  (-itA*^^)  $*m 

m  — 

v^  »  exp  (-itA**^^)  $*n 

The  behavior  for  t  of  these  three  functions  will  be  analyzed 
separately. 

Behavior  of  v^.  The  partial  wave  v^  has  the  representation 

(3.6)  Vjj^(t,x,y)  -  4i_(x,y,p,q)  exp(-i  t  a)^(p,q))  ll(p,q)  dpdq 

where 


(3.7) 


(Recall  X(p,q)  ■  w^Cp.q)  for  ±q  >  0.)  To  discover  the  behavior  of 
Vj^(t,x,y)  for  (x,y)  G  R*(d)  and  t  -►  “>  It  will  be  convenient  to  write 
<J)_(x,y,p,q)  In  a  way  that  puts  in  evidence  its  behavior  for  y  -*■  ±<».  To 
this  end  recall  that  by  [12,  (1.53)-(1.60) ] 

(3.8)  <j>_(x,y,p,q)  =  (27r)“*  c(«>)  (2q) e^*’  *  |pI.^) 

for  (p,q)  =  X^(p,X)  e  (and  hence  X  =  X(p,q)  -  c^(“)(|pl^  +  q^)). 

Moreover,  by  [12,  (4.3),  (4.6)  and  (8.41)]  one  can  write 


where 


(3.10)  T^(y,li,X)  -  T^(li,X)  4>^(y,y,X)  exp  {iy  q_(]j,X)} -♦•T^(U,X)  ,  y 
Similarly,  by  [12,  (4.1)  and  (4. 6)-(4.11) ]  one  can  write 


(3.11) 


f  00  ) 

'^+(y.y,^)  -  ®xp  {-iyq^(y,X)} 

+  R+(y»y»^)  exp  tiyq+(y»^)}] 


where 

' 

i^(y,y,^)  -  4>2(y*y*M  exp  {iyq^(y,X)}  -*■  i 
(3.12)-  •  y  -*■  +00. 

R^(y,y,^)  -  R^(y,X)  <l>i(y,y,X)  exp  {-iy  q_^(y,X)}->-R_j_(u,X) 


Combining  (3.6),  (3.8)  and  (3.11)  gives 
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Vjj^(t,x,y) 


c(oo)  p(<»)  (2Tr)"^/2 


exp  {i(x  •  p+y  q  -  tu)^(p,q)  }  I^(y,  |pj,X)  i!.(p,q)dpdq 


(3.13) 


+  c(«»)  P(~)*^^(2tt)  j  exp  {i(x  •  p- y  q  -  tu)^(p,q)  }  R.  (y,  jp  [ ,  X)  ll(p,q)dpdq . 

J  r 


It  is  natural  to  expect  that  in  R^(d)  the  partial  wave  v^(.t,K,y)  will 
propagate  as  t  -♦•  <*>  into  regions  where  y  is  large  and  hence  I_j_(y,  |p|,X) 
and  R^.(y,  1p1,X)  are  near  their  limiting  values.  Thus  the  representation 

(3.13)  suggests  the  conjecture  that 

(3.14)  ^  +v^(t,*)  inL2(R|(d)),  t  ■*  <=>> 

where  v^  and  v|  are  defined  by 


(3.15)  v^(t,x,y)  =  c(«>)p(<») (2tt)  j  exp  {i(x*p+y  q  -  tu)_|_(p,q)  })l(p,q)dpdq 


v^(t,x,y)  -  c(“)p(«)*''^  (2Tr)"^/*  I  exp  {i(x*p-y q  -  ta)^(p,q)}  x 


(3.16) 


R+(Ip|  >^)  Jl(p,q)dpdq 

-  c(“)p(<»)^^*  (2Tr)"^^2  I  {i(x*p+yq  -  tu)^(p,q)}  x 

"^+ 

^  R+(|p| .^)  A(p,-q)dpdq 


where  -C^  -  { (p,q)  :  (p,-q)  €  C^}  -  {(p,q)  :  q  <  -a|p|}.  Note  that  v® 
and  v^  are  waves  in  a  homogeneous  medium  with  density  p  (<»)  and  sound 
speed  c(®)  .  More  precisely, 
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(3.17) 


vS(t,  •) 


exp  (-it  c(“)  Ay^)  hf 


exp  (-i  t  c(“)  hi 


where  Ag  is  the  selfadjoint  realization  in  L2(R®)  of 

-A  -  -(3^/9xJ  +  3^/9x|  +  3^/3y^)  and  hj^  and  h^  are  the  functions  in 

L2(R^)  whose  Fourier  transforms  are 


(3.18) 


^^(p.q^  »  c(«>)  p(“)  iKp.q)  “  c(<»)  p(®) X+(P.q)  h_(p,q), 


1/2 


h£(p,q)  -  c(«>)  p(oo)  1/2  R^(jp  j,X)  A(p,-q) 


c(“)  p(®)  R^(|p|.A)  (1  -  X+(P.q))  h_(p,-q) 


Both  functions  are  in  Lj(R*)  because  h_  S  L2(R®)  and  |R^(|p,A)|  £p'^*('»)  by 
the  conservation  law  [12,  (1.46)].  Moreover,  supp  hj^  C  ~C^  and  hence 
the  theory  of  asymptotic  wave  functions  for  d'Alembert's  equation 
[9,  Ch.  2]  implies  that  v^(t,*)  -  0  in  L2(R^(d))  when  t  Combining 

this  with  (3.14)  gives 


(3.19)  Vj^(t,*)  v^(t,»)  in  L2(R|(d)),  t 

Now  consider  the  behavior  of  v^(t,x,y)  for  (x,y)  e  R^Cd),  t  -►  <». 
Combining  (3.6),  (3.9)  and  (3.10)  gives 


Vj^(t,x,y)  -  c(®)  p(<»)*^^  (211)”^^^  j  exp  {i(x*p+ yq_  -  tui^)  }  x 


(3.20) 


T^.(y,  |p  ,A)  «.(p,q)dpdq 


where  q_  -  q_(|p|,A),  -  u)^(p,q)  and  X  -  X(p,q)  -  u)^(p,q).  This 

representation  suggests  that 
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(3.21) 


Vn(t,»)  ~  Vo(t,*)  in  LjCR^Cd)),  t  “ 


where 


(3.22) 


v|(t,x,y)  -  c(“)  p(~)'/*  (2it)“®^^  j  exp  {i(x  •  p  +  yq_-  tto^)} 


r^(  1 P 1 » ^)  i!.(p,q)dpdq 


Now  the  mapping  (p,q)  -*•  (p,q')  =  X'(p,q)  (p.q_(  |p |  »t»J^(p»q) ) )  with 

domain  has  range  X' (C_j_)  =  =  rJ[(0),  Jacobian  3(p,q)/3(p,q ') 

m  c*(-<»)q'/c^(“)q  and  satisfies  (i)^(p,q)  =  a)_(p,q').  Thus  (3.22)  implies 


the  representation 


Vo(t,x,y)  ■  c(“)  p('»)^^^  (2ii)~®^^  f  exp{i(x*  p  +  yq'  -tu  )}  x 


(3.23) 


^  T^(|p1  .“l(p.q'))  X'(p.q)  (c^(-“)q7c*(“)q)dpdq' 


where  q  *  q(lpl,q*)  *  /a* (  p  ^  +  q’^)  +  q'^ .  Note  that 


(3.24) 


v|(t,*)  -  exp  (-itc(-'»)  aJ/^)  h| 


where  h|  e  L2(R^)  has  Fourier  transform 


(3.25)  hj(p,q')  ■  c(«)  p(®)  T^(  jp  ]  (p,q'))  X.(p,q)  (c^(-®)q7c^(«>)q)  . 

Since  supp  h^  C  the  results  of  [9,  Ch.  2]  imply  that  v|(t,  •)  '•  0  in 
L2(R*(d))  when  t  -►  <*>.  Combining  this  with  (3.21)  gives 


(3.26) 


Vj(t,*)  ~  0  in  L2(R®(d)),  !-►<». 


Analogous  conjectures  concerning  Vjjj(t,*)  and  Vj^(t,*)  will  now  be 
formulated.  Only  the  main  steps  of  the  calculations  will  be  given  since 
the  method  is  the  same  as  for  V0(t,*). 
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Behavior  of  v  .  v  has  the  representation 
_ m  ni 

(3.27)  v^(t,x,y)  =  <P_(x,y,p,q)  exp  (-1 1  OJj^Cp.q))  m(p,q)dpdq, 

•'  Co 

by  (3.5),  where 

(3.28)  <PJx,y,p,q)  -  (2TT)-^c(<»)(2q)^/^ 

for  (p,q)  ■  X(|(p,X)  €  C,  (and  hence  X  ■  X(p,q)  *  w^(p,q)).  Moreover 
(see  [12,  (A.ia)-(4.2A)]), 


W/2 

(3.29)  'Po(y.U.^)  “  [47r<^(;j,X)J  ^^(y.vi.X)  exp  (yqUy,X)) 


where 


(3.30)  To(y,y,X)  -  To(y,X)  (|»3(y,y,X)  exp  (-yqJ^(y,X))  -►  T5(y,X),  y  ->■  -« 


and 


(3.31) 


where 


'j^o(y»^^.^)  “  [4u<^^,X)]  ^  {Io(y»M.A)  exp  {-iy  q^(y,X)} 
+  Ro(y,y,X)  exp  {iyq^(y,X)}] 


(3.32)  \ 


I,  (y,y,X)  -  4i2(y.M.^)  exp  tly  q^,(y,X)}  -*■  1 


Ro(y*W*^)  -  <(ij(y,y,X)  exp  {-iy  q^(y,X)}-^Rj(y,X) 
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Combining  (3.27),  (3.28)  and  (3.29)  gives 


v^(t,x,y)  =  c(«>)  p(“>)*^^  (2ii) 


exp  {l(x  •  p  -  tco_j_)  }  X 


(3.33) 


X  TQ(y,|p|,X)  exp  (y q^)  m(p,q)dpdq 


where  u)^  -  aj^(p,q),  X  *  w^(Pfq)  and  q^  »  92(|p|»^)*  Since 


(3.34) 


T.(y,  1p1  ,X)  exp  (y  ql(|pl  ,X))  ->-0,  y 


equation  (3.33)  suggests  that 


(3.35) 


v^(t,*)  ~  0  in  L2(R^(d)),  t  -*•  «. 


Similarly,  combining  (3.27),  (3.28)  and  (3.31)  gives 


v^(t.x,y) 


c(oo)  p(oo)*/2  (2ir) 


1/2 


exp  {i(x*p+yq-  taj_|_) }  x 


X  l5(y, |p|,X)  m(p,q)dpdq 


(3.36) 


+  c(«’)  p(®)  (2ir)  I  exp  {i(x  •  p  -  yq  -  tu^)} 


X  Ro(y,|p  ,X)  m(p,q)dpdq 


which  suggests  that 


(3.37)  v^(t,-)  ~  v«(t,*)  +  vi(t,*)  in  L2(R^(d)),  t  « 


where  v®  and  vj;  are  defined  by 
mm 


(3.38) 


v^(t,0  -  exp  (-itc(®)  aJ/*)  hjjj 


^i(t»*)  “  ®*P  (-itc(«)  Aj/*)  h' 


> 
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and  h  and  are  the  functions  In  L2(R^)  whose  Fourier  transforms  are 
h^(p»q)  *  c(<»)  m(p,q)  -  cC")  p(«>) Xo  h_(P»<l). 

(3.39)  ' 

h^(p.q)  ■  c(®)  P(“)  RjCIpI  ,X)  ra(p,-q) 

-  c(«»)  p(“)  Ro(|p1  ,A)  (1  -  Xo(P*q))  h_(p,-q)  . 

Note  that  supp  h^  C  and  hence  v*(t,*)  ~  0  In  L, (R?(d))  when  t 
Combining  this  with  (3.37)  gives 

(3.40)  In  L2(R^(d)),  t 

Behavior  of  v  .  v  has  the  representation 
_  n  n 

(3.41)  V  (t,x,y)  »  4i^(x,y,p,q)  exp  (-1 1  aj_(p,q))  n(p,q)dpdq, 

7c_ 

by  (3.5),  where 

(3.42)  <|)_(x.y.p,q)  -  (2tt)’‘c(-<»)  (2  jq  ] )  ^^^  e^P*’'  ip.(y.|pl,A) 

for  (p,q)  *  X_(p,X)  €  C_  (and  hence  X  »  X(p,q)  ■  u)f(p,q)  .  Moreover 
(see  [12,  (4.5)  and  (4.12)-(4.17) ]) , 

(3.43)  '|J_(y,P,X)  -  [4^q^(^|'X)]  T_(y,u,X)  exp  {lyq^(y,X)} 

where 

(3.44)  T_(y,y,X)  -  T_(y,X)  (y,u,X)  exp  {-ly  q^(y ,X) }  ■>  T_(u,X) ,  y -*•  « 
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and 


(3.45) 


’J^_(y.U.X)  »  [4Vq^^X)]  ^  II_(y.U.X)  exp  {iy  q_(y,X)} 
+  R_(y,U,X)  exp  {-iy  q_(M,A)}] 


where 


I_(y.U,A)  «  <j)3<y,U,A)  exp  {-iyq_(p,X)}  -*•  1 


(3.46)  ( 


R_(y.y.A)  -  R_(y,X)  <J>4(y,y,X)  exp  {iy  q_(y,X)}-^R_(y,X) 


Combining  (3.41),  (3.42)  and  (3.45)  gives,  after  simplification  using 
q_(|p|»“!(p.q))  “  (q^)^^“  -q  for  (p,q)  e  c_. 


Vjj(t,x,y) 


(3.47) 


c(-“)  p(-“>)  (Itt)  exp  {i(x»p+yq-  tu)  ) }  x 

•'C 


^  f_(y»|P  .A)  n(p,q)dpdq 


+  c(-«»)  p(-«>)*^^  (2it)~®^^  exp  {i(x*  p-yq  -  tu)  )} 

■'C_ 

R_(y»  1p1  .A)  n(p,q)dpdq. 


This  suggests  the  asymptotic  behavior 

(3.48)  L2(R®(d)),  t  -►  », 


where 


(3.49) 


v®(t,*)  -  exp  (-itc(-«»)  A^*)  h^ 
v^(t,*)  -  exp  (-ltc(-«»)  aJ^*)  h^ 
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and  h  and  are  the  functions  whose  Fourier  transforms  are 
n  n 


(3.50) 


h^(p.q)  =  c(-“)  p(-“) n(p.q)  =  c(-»)  p(-“)  X_(p.q)  h_(p,q), 


I  h^(p,q)  =  c(-“)  p(-“)  R_(1p1  .X)  n(p,-q) 

=  c(-<»)  p(-oo)  1/2  R_(  jp|  ,X)  (1  -  x_(P.q))  h_(p,-q). 

Note  that  supp  h^  C  and  hence  v^(t,*)  ~  0  in  when  t 

Combining  this  with  (3.48)  gives 

(3.51)  inL2(Rl(d)),  t  - 

Finally,  combining  (3.41),  (3.42)  and  (3.43)  gives 


v^(t,x,y) 


=  c(-<»)  p(-“)*^^  (2Tr)“®^^ 


exp  {i(x  •  p  -  y  q_^  -  tu)_)  }  x 


(3.52) 


^  T^(y, 1p|,^)  n(p,q)dpdq 

where  u)_  ■  u_(p,q),  X  ®  u)^(p,q)  and  q_j_  »  q^(|p|»X).  This  suggests  that 

(3.53)  ^n^*^**^  ~  v2(t,*)  inL2(R®(d)),  t  -><» 


where 


v^(t,x,y) 

(3.54) 


c(-oo)  p(-<»)  (2it)""*^^ 


exp  {i(x  •  p  -  y  q^.  -  tii)_ 


)} 


T_(  IpI  ,X)  n(p,q)dpdq. 


Now  the  mapping  (p,q)  -*■  (p,q')  -  X"(p,q)  •  (p,-q^(p,X(p,q)))  maps  C_ 
onto  X"(C_)  ■  -C^,  has  Jacobian  3(p,q)/3(p,q ’)  ■  c^('“)q '/c*(-“>)q  and 
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satisfies  u_(p,q)  *U)^(Pfq').  Thus 


(3.55) 


v^(t,*)  *  exp  (-itc(»)  aJ^^) 


where  h^  has  the  Fourier  transform 
n 


h^Cp.q')  “  c(-“)  p(-oo)  1-^2  T_(|p|,u2(p,q'))  n(p,q(p,q')) 


(3.56) 


X  (c^(“)q7c^(-<»)q) 


Moreover,  supp  C  and  hence  v^(t,*)  ~  0  in  L2(R^(d)),  t 
Combining  this  with  (3.53)  gives 


(3.57) 


v^(t,*)  ~  0  in  L2(R^(d)),  t 


The  asymptotic  behavior  of  Vj(t, •)  for  t  “  may  be  obtained  from 
the  three  cases  analyzed  above  by  superposition,  equation  (3.4).  Thus 
equations  (3.19),  (3.26),  (3.35),  (3.40),  (3.51)  and  (3.57)  imply 


v^(c,*)  +  v^(t,»)  in  L2(R{(d)) 


(3.58) 


Vj(t,*) 


v“(t.-) 


in  L2(R!(d)) 


t 


On  combining  this  with  the  definitions  of  v^,  v^  and  v^,  equations  (3.17), 
(3.18),  (3.38),  (3.39),  (3.49)  and  (3.50),  one  is  led  to  formulate 
Theorem  3.1.  For  every  h  €  JC  let  Vj(t,»)  be  defined  by 


(3.59)  v®(t,x,y) 


fexp  (-itc(®)  aJ^^)  h^(x,y),  (x,y)  e  R^(d) , 


exp  (-itc(-«)  Ai'^*)  h“(x,y),  (x,y)  e  Rf  (d) , 


where  h'*'  and  h~  are  the  functions  in  L2(R^)  whose  Fourier  transforms  are 


given  by 


1 
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(3.60) 


h'''(p,q) 


c(«>)  P(") h_(p,q)  ,  (p,q)  €  R_J., 


,  (p,q)  e  R^, 


and 


(3.61) 


vs  . 

h  (p.q) 


0  ,  (p,q)  e  R^, 

c(-“)  p(-“)  h_(p,q) ,  (p,q)  e  R^. 


Then 


(3.62) 


llm  1  v^(t,  •)  -  v°(t,  on 


K 


0. 


Theorem  3.1  implies  corresponding  asymptotic  estimates  for  the 
free  component  u^(t,*)  =  P^u(t,*)  =  Re  {v^(t,*)}  of  the  acoustic  potential 
u(t,*).  Indeed,  if  u^(t,*)  is  defined  by 


(3.63)  u“(t,*)  a  Re  {vj(t,‘)} 

then  Theorem  3.1  and  the  elementary  inequality  |Re  z|  <  |z|  imply 
Corollary  3.2.  For  all  h  €  3C  one  has 


(3.64) 


llm  9  U£(t,  •) 

f*co 


U£(t,  Olljf  =  0. 


If  the  initial  state  h  has  derivatives  in  Jf  then  Uj(t,*)  and 
u®(t,  •)  have  the  same  derivatives  in  and  (3.64)  can  be  strengthened  to 
include  these  derivatives.  In  particular^  one  has 

Corollary  3.3.  For  all  h  €  L^(R^)  -  D(A^/2)  one  has 


(3.65) 


llm  JD.  u^(t,*)  -  D.  ul(t,*)lljf  “  0,  j  =  0,1, 2, 3, 
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where  *•  9/3t,  “  3/3xi,  *  3/3x2  “  3/3y.  Equation  (3.65) 

is  equivalent  to  convergence  in  energy: 

(3.66)  lim  E(uj  -  U£,R^,t)  =  0. 

t-xn 

Corollary  3.3  can  be  proved  by  applying  the  method  of  this 
section  to  the  derivatives  v^(t,*)  (j  *  0,1, 2, 3)  which  are  given  by 

integrals  of  the  same  form  as  (3.2).  Detailed  proofs  for  the  case  of  the 
Fekeris  profile  were  given  in  [10]. 

Proof  of  Theorem  3.1.  The  remainder  of  this  section  is  devoted 
Co  the  proof  of  Theorem  3.1.  The  decomposition  (3.3)  is  used  for  the 
proof.  Moreover,  for  brevity,  only  the  asymptotic  equality  (3.53)  for 
v^(t,*)  is  proved.  The  remaining  five  cases,  namely  (3.1A),  (3.21), 
(3.35),  (3.40)  and  (3.48)  can  be  proved  by  the  method  used  for  (3.53). 

As  a  first  step,  (3.53)  will  be  proved  for  the  special  case  of  n(p,q) 

€  Cq(C_),  the  set  of  continuous  function  with  compact  supports  in  the 
open  cone  C_.  The  general  case  will  then  be  proved  by  using  the  fact 
that  Cq(C_)  Is  dense  in  L2(C_). 

For  functions  n(p,q)  €  Cg(C_)  the  integrals  defining  v^  and 
converge  point-wise,  as  well  as  in  7C,  and  one  can  write 


where 
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(3.68)  w(y,p,t) 


'0 

exp  {-i(yq_^  +  ta)_)}  lT_(y,  |p|  .tuf)  -  T_(  |p|,a)f)  Jn(p,q)dq 

J  .00 


Parseval's  formula  in  L2(R^),  applied  to  (3.67),  gives 


(3.69) 


Vj^(t,x,y)  -  v^(t,x,y)  (^dx 


c^(-«)  p(-»)(27r)~‘[  |w(y,p,t)  1^  dp. 


On  integrating  this  over  y  >  d  one  finds 

The  last  relation  implies  that  to  prove  (3.53)  it  is  sufficient 
to  prove  that 

(3.71)  w(-,t)  ^  0  in  L2(R|(d)),  t  ^  ». 

To  this  end  it  will  be  convenient  to  change  the  variable  of  integration 
in  (3.68)  from  q  to  U)  *  w_(p,q)  »  c(-<»)  /|p  |  ^  +  9^  •  Solving  this 
equation  for  q  <  0  gives  q  ■  -(u)^  c  ^(-”)  -  lp|^)*^^  *  “9_(|p|>w^)  with 
0)  >  c(-“)|p|.  Hence  (3.68)  can  be  written 


(3.72) 


where 


w(y,p,t) 


IpI 


exp  (-itu))  W(y,p,u))du) 


W(y,p,w) 

(3.73) 


-2 

C 


(-*)  exp  {-iy  q^(  lpl,u)^)}[T_(y,  jp|  ,u)^)  -  T_(  jp|  ,u)^)  ]n(p,-q_(  jp]  ,a)^))a) 


q_(|p|  ,0)^) 
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i 


I 


I 


♦ 


The  assumption  that  n  €  Cg(C_),  together  with  (3.44),  Implies  that 
W  e  Cg(R  X  r)  where  T  »  {(p,w)  :  (n  >  c(-«)|p|}.  Moreover,  by  a  standard 
partition  of  unity  argument,  one  may  assume  without  loss  of  generality 
that 

(3.74)  supp  W(y,*)  C  {(p,t»»)  :  jpj  <  Po  and  0  <  Uj  <  oj  <  ajj} 
for  all  y  €  R  where  (jj^j  >  c(-<»)pg.  This  in  turn  Implies  that 

(3.75)  w(y,p,t)  -  I  *  exp  (-itu))  W(y,p,u))da) 

and 

(3.76)  supp  w(y,*,t)  C  B(po)  -  {p  :  [p]  <  p^^} 
for  all  y  €  R  and  t  €  R.  Thus 


(3.77) 


lw(«,t)|*  ■  I  f  |w(y,p,t)|2  dpdy 

-  [  [  Iw(y,p,t) I*  dpdy 

^B(pJ 


+  j  f  lw(y,p,t)l^  dpdy 

^B(p,) 


for  any  y^  >  d.  The  proof  of  (3.71)  will  be  derived  from  (3.77)  and  the 
following  two  lemmas. 

Lemma  3.4.  Let  n€  Cg(C_)  and  assume  that  (3.74)  holds.  Then 
for  each  d  €  R,  y^  >  d  and  Pg  >  0  one  has 
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(3.78)  llm  w(y,p,t)  -  0, 

t-«s 

uniformly  for  all  (y,p)  S  [d.yj]  x  BCp^) . 

Lemma  3.5.  Under  Che  hypotheses  of  Lemma  3.4  there  Is  a  constant 
C  *  C(n)  such  that 

(3.79)  |w(y,p,t)i  <  C  y'““ 

for  all  y  >  0,  p  e  B(pfl)  and  t  e  R,  where  a  >  3/2  is  the  constant  of 
condition  (1.6). 

Proof  of  Lemma  3.4.  The  proof  is  based  on  a  well-known  proof  of 
the  Riemann-Lebesgue  lemma.  Note  that  by  (3.75)  one  has 


w(y,p,t) 


(3.80) 


■(i>j-(tr/t) 

exp  (-iu  t)  W(y,p,aH-(TT/t))da) 

''a)j-(Tr/t) 


exp  (-i  t  u)  [W(y,p,(D)  -  W(y,p,u)+(u/t))  ]du) 


1  f“o 

-5  exp  (-itw)  W(y,p,a)+(ir/t))dw 

•'a)o-(TT/t) 


+  4  I  ^  exp  (-itu)  W(y,p,urf(iT/t))dio. 

■'a)j-(Tr/t) 

The  limit  relation  (3.78)  is  obvious  from  (3.80)  and  the  continuity  of 
W.  The  uniformity  of  the  limit  follows  from  (3.80)  and  the  uniform 
continuity  of  W  on  compact  subsets  of  R  x  r. 

Proof  of  Lemma  3.5.  Note  that  by  (3.72),  (3.73)  and  (3.44)  one 


has  the  estimate 


I 


9 


9 
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(3.81) 


|w(y,p,t) 


(y.  IpI.w^)  exp  {-iy  «4(1p1.w^))-  l|  M(p,u))du 


for  all  y,  p  and  t  where 


(3.82) 


M(p,a3)  =  c“®(-")  1t_(  IpI.oj^)  I 


n(p,-q_(  |pl,w*))a) 


<!_(  IpI.w*) 


is  continuous  for  1p1  f  Po»  <  (o  <  Wj.  It  follows  that  for  any  y  €  R, 
p  €  B(pg)  and  t  €  R 

(3.83)  lw(y,p,t)l  <  Mg  sup  l<l>i  (y,  jp  |  ,<»)^)  exp  {-iy  q^(  Ip  |  ,U)^)  }  -  l| 

where  Mg  ■  Mg(n)  *  sup  M(p,a))  and  the  suprema  are  taken  over  all  |p|  <  Pg 
a".d  Ug  <  03  <  Uj.  The  proof  of  (3.79)  will  be  based  on  (3.83),  conditions 
(1.3),  (1.4),  (1.6)  on  p(y)  and  c(y)  and  the  proof  in  [12]  of  Theorem 
2.1;  see  [12,  p.  27ff].  Note  that  because  of  the  continuity  of  w  it  will 
suffice  to  prove  (3.79)  for  all  y  >  y^  where  yj  ■  yi(n)  is  a  positive 
constant. 

The  solution  <|»i(y,M,A)  with  X  >  c^(-<»)ij^  >  c^(<»)y^  satisfies 
[12,  p.  27ff] 

(3.84)  <l»i(y,M,X)  -  exp  {lyq^(y,X)}  (n^  +  Hj) 

where  n  -  (rij,n2)  Is  characterized  on  y  >  yj  as  the  unique  solution  of 
the  Integral  equation  [12,  (2.64)]  which  can  be  written 


(3.85) 


n  -  n®  +  K(w,X)n,  n®  -  (i,o). 


The  kernel  K(y,y',u,X)  -  (K^j (y,y' ,u,X))  is  defined  by 


I 
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(3.86) 


and 


Kjj(y,y',U,A) 


0  ,  Vity'  y« 


y'  >  y. 


(3.87)  K2j(y,y’.P.A)  = 


0  .  yi  <  y'  <  y. 

|-exp  {-2i(y-y’)q^(U,X)  E2j(y',p,X),  y'  >  y. 


where  [12,  p.  27] 


(3.88) 


E(y,y,X)  -  B”*(y,X)  N(y,y,X)  B(y,X). 


From  these  relations  one  has  [12,  (2.64)] 


'J>i(yty.^)  exp  {-iyq^(y,X)}  -  1  »  m  +  -  1 


(3.89) 


fCD  000 

-J  Ejj(y',y,X)  nj(y')dy’  -  exp  {-2i(y-y')q^(y,X)}  (y ' ,M,X)nj (y ')dy ' 


and  hence 


2  f®® 

(3.90)  l<(ii(y,y,X)  exp  {-iy  q^(y,X)}- l]  <  f  j  lE.j^(y' ,y,X)  ]  Irijj(y')  jdy ' . 

j,k-l  Jy 

Using  [12,  (2.66)]  and  the  continuity  of  B(p,X)  on  X  >  c^(-«')y^  it  can 
be  shown  that 

(3.91)  lK(y,X)l  <  1/2  for  0  <  y  <  Pj ,  <  X  < 

provided  yj  ■  yi(n)  is  large  enough.  Thus 
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(3.92) 


Ini  < 


-  1  -  IKI  1  -  IKI  - 


<  2 


for  0  <  y  <  Pq,  ojg  <  X  <  uif.  Combining  this  with  (3.90)  and  [12,  (2.66)] 
gives 


(3.93) 


I 

^.(y.y.X)  exp  {-iy  q.(y,X)}- l|  <2  \  |E.(y',y,A)|  dy' 

h  j.k-x  ^ 

<  C,  I  Ip(y')  -  p(“) 1  dy' 

*  V 


for  the  same  values  of  y  and  X  where  Cj  and  depends  only  on  Pj,  oj^ 
and  (Uj.  It  follows  from  (3.93)  and  (1.6)  that 

(3.94)  sup  |4'i(y,  1p1.<jJ^)  exp  {-iyq^dpl.w^)}  -  l|  <  Cj  y*"® 

for  all  y  >  yj(n)  where  C,  depends  only  on  Pj,  u)^  and  ui^  (i.e.,  n)  . 
Combining  (3.94)  with  (3.81)  gives  (3.75). 

Proof  of  Theorem  3.1  (completed) .  Lemma  3.5  Implies  that 

(3.95)  r  [  |w(y,p,t)|*  dpdy  <  (Trp*  C^/2a  -  3)  y, 

Vo  ^B(Po) 

where  3  -  2a  <  0.  Thus  given  any  e  >  0  there  is  a  y^  ■  yg(£,n)  such  that 


(3.96) 


B(P,) 


lw(y,p,t) 1* 


dpdy  <  e 


for  all  t  €  R.  Equation  (3.77),  together  with  Lemma  3.4  and  the 
estimate  (3.96),  Imply  that 
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(3.97) 


11m  sup  l|w(*,t)ll‘  <  e. 


L2(R^(d))  - 


Since  e  >  0  Is  arbitrary  this  Implies  (3.71),  as  required. 


The  arguments  given  above,  applied  to  Vj^^,  v^  and  v^^,  show  that 
the  conclusions  of  Theorem  3.1  hold  for  all  h  such  that 
h_  €  Cj(C^  U  Cq  U  C_) .  Moreover,  this  set  Is  dense  In  L2(R*)  and  hence 
U  Cj  u  CJ  Is  dense  In  =>  X  by  [12,  Cor.  8.12].  These 
facts  can  be  used  to  extend  (3.62)  to  all  h  e  X  because  the  mappings 
U(t)  ;  X^  -►  L2(R®)  and  Uj(t)  :  X^  ->■  L^iR^)  defined  by  U(t)h^ 

■  exp  (-ItA*^^)  hj  and  Uj(t)hj^  =  V£(t,»)  are  uniformly  bounded  for  all 
t  €  R  (see  [11,  p.  32]).  The  density  argument  needed  to  extend  (3.62) 
to  all  h  €  JC  has  been  given  In  many  places;  see,  for  example,  [9,  Ch.  2] 
or  [10,  p.  260]. 
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§4.  Transient  Guided  Waves. 

The  asjnnptotic  behavior  for  t  “  of  the  guided  component 
Ug(t, •)  »  Pg  u(t, •)  is  derived  in  this  section.  Ug(t,*)  is  a  sum  in  3f 
of  mutually  orthogonal  partial  waves  Uj^(t,*)  ®  Pj^  u(t,*)  *■  Re  {vj^(t,*)}, 
1  5  k  ^  Nq.  The  starting  point  for  the  analysis  is  the  Integral 
representation 


(4.1) 


Vj^(t,x,y)  »|  ip]^(x,y,p)  exp  {-it  Wj^(IpI)  hj^(p)dp 


where 


(4.2) 


\(p)  =  I  ^  '('j^(x,y,p)  h(x,y)  c"*(y)  p"^(y)dxdy 


and  the  Integrals  converge  in  3C  and  L2((ij^)»  respectively.  The  Integral 
in  (4.1)  can  be  written 


(4.3)  Vj^(t,x,y)  ■•^1  exp  {i(x  •  p -  t  ujj^(  [p  | ) ) }  (//^^(y.p)  ^(p)dp. 


This  is  an  oscillatory  Integral  that  can  be  estimated  by  the  method  of 
stationary  phase  when  hj^  €  C*(fij^)  .  To  apply  the  method  define 


r  »  /t*  +  Ixl 


(4.4) 


t  -  r  Co.  *1  ■  r  Cl,  Xj  -  r  C2 


C  -  (Co.Ci.Cj)  €  S*  C  R» 
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where  denotes  the  unit  sphere  In  R^.  Then  (4.3)  takes  the  form 


(4.5) 


Vj^(t,x,y) 


exp  {irej^(p,C)}  gj^(p,y)  dp 


where 


(4.6) 


0jj(p.O  -  CiPi  +  CaP2  -  Co  ‘^kOp|) 
glt(p.y)  “  (2ir)~*  t(»|^(y,p)  hj^(p) 


Estimates  for  large  r  of  the  Integral  In  (4.5)  are  needed  that  are 
uniform  for  (C,y)  In  compact  subsets  of  x  R.  Such  estimates  are 
provided  by  a  version  of  the  method  of  stationary  phase  due  to 
M.  Katsumura  [5].  A  form  of  Matsumura's  results  applicable  to  (4.5)  ras 
presented  In  [11,  Appendix].  This  result  Is  applied  below  to  estimating 

The  phase  function  has  a  point  of  stationary  phase  If  and  only 

if 

(4.7)  Ujj(Ip1)p/  IpI  “  x/t 


where 


(4.8)  \(!p|)  -  u)^(1p1) 

Is  the  group  speed  associated  with  the  dispersion  relation  o)  -  b>j^(|p|)> 
It  will  be  assumed,  for  brevity,  that  l^]^(Vi)  Is  a  monotone  decreasing 
function  that  maps  0^^  onto  (c^,c (<»)).  In  this  case  (4.7)  has  a  unique 
solution  If  |x|/t  lies  In  the  range  of  U|^(|p|);  that  Is 
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(4.9) 


<  |x|/t  <  c(®), 


and  no  solution  otherwise.  The  solution  is  given  by 


(4.10) 


P  ”  Qv( l*l/t)  x/ |x 


where  is  the  inverse  function  to  Uj^.  By  calculating  the  Hessian  0^ 


one  can  show  that 


(4.11) 


Idet  9I'(p,C)|  *  t^  U.(|p|)  |U'(|p|)  |/|p| 


and  sgn  0j^(p,5)  ■“  0.  In  particular,  each  point  of  stationary  phase  is 
non-degenerate  and  makes  a  contribution 


(4.12)  V.  (t,x,y,p)  =■ 


(p|^''^exp  {i(|x(  ip|-ta)^.(|p|))}  i//^(y,|p|)h^^(p) 
t{Uv(lpl)lu'(lpl)l}'/2 


to  the  integral  in  (4.3),  where  p  is  given  by  (4.10).  For  jxj/t  outside 
the  interval  (4.9)  there  is  no  point  of  stationary  phase.  Thus  the 
stationary  phase  approximation  to  Vj^(t,x,y)  is  given  by 


(4.13) 


vr(t,x,y)  “  x(|x|/t)  v"(t,x,y,Q^( jxj/t)  x/jxl) 


where  X  the  characteristic  function  of  the  interval  (c^,c(“>))  and  one 


Theorem  4.1.  For  all  h  e  5C  such  that  hj^  e  Cg(fij^)  there  exists  a 
constant  C  •  Cj^(h)  such  that 


(4.14) 


|vj^(t,x,y)  -  v^(t,x,y)|  <  C/t' 


for  all  t  >  0,  X  €  -  {O}  and  y  e  R. 


Theorem  4.1  can  be  proved  by  application  of  Theorems  A.l  and  A. 2 


of  [11].  The  proof  that  has  the  required  p-derlvatlves  Is 

lengthy  but  straightforward  and  will  not  be  given  here.  If  hj^  is  not  a 
smooth  function  then  Theorems  A.l  and  A. 2  are  not  applicable  and  the 
estimate  (4.14)  may  fail.  However,  the  definitions  (4.12)  and  (4.13)  are 
meaningful  for  all  h  €  K  and  one  has 

Theorem  4.2.  For  all  h  e  JC,  all  t  >  0  and  k  =  l,  2,  3,***  one  has 

(4.15)  v^(t, 

and 


(4.16)  * 

CO 

Moreover,  the  mapping  t  Vj^(t,*)  is  continuous  from  to  3C  and 

(4.17)  11m  lvj^(t,*)  -  v^(t,*)*jj,  -  0. 

The  proofs  of  these  properties  are  the  same  as  those  for  the 
Pekerls  profile,  given  in  [10],  and  are  not  reproduced  here.  On 
defining 

(A. 18)  u^(t,x,y)  -  Re  {v^(t,x,y)} 


one  also  has 

Corollary  4.3.  For  all  h  €  JC  and  k»l,  2,  3,*»*, 


(4.19) 


11m  luj^(t,*)  -  u^(t,»)l3f 

t-*«9 


0. 


If  h  €  L2(R*)  then  u^(t,*)  €  L2(R*)  and  asymptotic  wave  functions 
for  the  first  derivatives  of  u^  can  be  constructed.  Indeed,  if 
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e  then  the  first  derivatives  of  Vj^  are  given  by 


°t 

(4.20) 

Dj  Vj^(t,x,y) 
Dy  vj^(t.x.y) 


x 

2tt 


27r 


iTT 


exp  {i(x  •  p  -  t  aj^.(lp|))}(-ia)j^(  |p|))ipj^(y,p)hj^(p)dp, 


exp{i(x*p-t  oJj^CIpI))  (iPj)  \(p)dp  (j  = 


exp  {i(x  *  p  ~  t  aJj^dpj))}  Dy  iJ^j^Cy.p)  hj^(p)dp. 


1,2), 


These  integrals  have  the  same  form  as  the  integral  (4.3)  for  v^^.  The 
corresponding  asymptotic  wave  functions  are  defined  by 


(4.21) 


v^o(t,x,y,p)  -  (-i  Wj^dp)))  v^(t,x,y,p). 


v*j(t,x,y,p)  =  (ipj)  v^(t,x,y,p)  (j  =  1,2), 


v^3(t,x,y,p)  =  Dy  v"(t,x,y,p) ,  and 

v“j(t,x,y)  =  xdxi/t)  v“j(t,x,y,Q^(  ix|/t)  x/jxj) 

for  j  ■  0,1, 2, 3.  The  analogue  of  Theorem  4.2  is 

Theorem  4.4.  For  all  h  e  lJ(r’),  all  t  >  0  and  k  »*  1,2,3, •••  one 

has 


(4.22) 


Vy(t,')  e  LjCR’),  j  -  0,1, 2. 3, 
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and 


(4.24) 


llm  DD 

t-wo 


(R') 


0,  j  =  0,1. 2, 3. 


The  proof  of  Theorem  4.4  is  the  same  as  for  the  special  case  of 
the  Pekerls  profile  which  was  treated  In  detail  In  [10]. 

The  preceding  discussion  was  restricted  to  the  special  case 
where  Uj^(lJ)  is  monotonic.  If  has  a  finite  number  of  maxima  and 

minima  there  are  a  corresponding  number  of  points  of  stationary  phase 
and  the  form  of  the  asymptotic  wave  function  is  more  complicated  but 
still  tractable.  In  the  case  of  the  Pekerls  profile,  treated  in  [10], 
there  are  two  points  of  stationary  phase.  Cases  that  lead  to 
infinitely  many  stationary  points  have  not  yet  been  encountered.  They 
would  require  additional  analysis. 


§5.  Asymptotic  Distributions  of  Energy  for  Large  Times. 

The  total  energy  of  the  acoustic  field  u(t,x,y),  given  by  (2.14), 
is  constant  for  t  >  T.  The  same  is  true  of  the  partial  waves  u^,  u^  and 
Uj^,  k  -  1,  2,  Moreover,  it  was  shown  in  [12]  that  .Pj,  *  •*) 

is  a  complete  family  of  orthogonal  projections  in  K  that  reduces  A.  It 
follows  that 


No-1 

(5.1)  lA^'^  hl^  -  lA*/^  h^O J  +  I  lA*/^ 

k“l 


which  may  be  Interpreted  as  an  energy  partition  theorem.  The  partial 
energies 


(5.2) 


E(uj,RSt)  -  lA^^  hjSS 

E(u^^,R^t)  -  lA^^^hj^il^,  k  -  1,  2.  3,*“,  Nq-1, 


can  be  calculated  from  the  source  function  f(t,x,y)  and  the  normal  mode 
functions.  The  relationship  between  h  and  f  is  given  by  (2.7),  which 
implies  that 


(5.3)  h_(p,q)  -  i  (p,q)  f_(X*^^  (p,q).P.q) 


where 


(5.4)  f_(w,p,q) 


exp  (IwT)  (p_(x,y,p,q)  f (r.x.y)  c“^ (y)  p“' (y)dydxdT 


and  similarly 
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(5.5) 


where 


\(p)  ■  i  \*(|p!)  fj^c^dph.p) 


(5.6)  f^(w,p)  -  ll  ^  exp  (iuT)  (j>j^(x,y,p)  f(T,x, 


y)  c  *(y)  p“'(y)dydxdT. 


It  follows  chat  the  partial  energies  are  given  by 


(5.7)  E(u.,R*,t)  -  [  If  (X^^*  (p,q)  ,p,q)  I*  dpdq 

*  Jb*  “ 


and 


(5.8) 


E(Uj^,R*,t) 


I  |fi.(“i,(|p|).p)  I*  dp.  k  -  1.  2, 


for  every  t  >  T. 

The  theorems  of  §3  and  §4  make  it  possible  to  calculate  asymptotic 
distributions  of  energy  in  bounded  and  unbounded  subsets  of  R^.  Only  the 
principal  results  are  formulated  here.  The  proofs  are  omitted  since 
they  are  the  same  as  chose  for  the  Pekerls  profile  which  were  given  in  [10] . 

The  notation 


(5.9)  E*(u,K)  «  lim  E(u,K, t) 

will  be  used  whenever  the  limit  exists.  A  first  result  is  the  transiency 
of  all  waves  with  finite  energy  in  stratified  fluids: 

(5.10)  E*(u»K)  ■  0  for  all  compact  sets  K  C  R®. 
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The  Free  Component.  The  results  of  §3  imply  that  in  each  half- 
space  R^(d)  (resp.  R^Cd))  behaves  like  a  wave  in  a  homogeneous 

medium.  Thus  if  (resp.  C  )  denotes  a  cone  in  R^(d)  (resp.  R2(d)) 
then  Corollary  3.3  and  the  results  of  [9]  imply 


(5.11) 


^(u,,C^)  -  c^(±“)  p(±®)  [  (IpI^  +  q*)  |h_(p,q)l^  dpdq 

Jc- 

-  P(±“)  j  ^  lf_(u^(p,q) .p.q) I*  dpdq. 


It  follows  that  if 

(5.12)  S  -  {(x.y)  :  dj  <  y  <  dj} 
is  a  slab  then 

(5.13)  E“(uj,S)  =  0. 

The  Guided  Component.  Consider  the  family  of  cones  defined  by 

(5.14)  C(E,d)  -  {(x,y)  ;  ly  -  d|  <  e|x|} 
where  e  >  0.  Then,  in  contrast  to  (5.11),  one  has 

(5.15)  E”(uj^,C(e,d))  -  E"(uj^,R®)  -  hj^l  J 

for  every  e  >  0  and  k  -  1,  2,  3,*»*  (see  [10,  Theorem  5.5]).  Finally, 
if  S  is  the  slab  defined  by  (5.12)  then  one  can  show  by  means  of  Theorem 
4.4  that 


E“(u^^tS)-|^  lfk(<\(|p|).p)  I' 


j  'l'k(y.P)  c"*(y)  p"‘(y)dy 


dp. 


(5.16) 
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§6.  Seml-Inflnlte  and  Finite  Layers. 

The  preceding  analysis  Is  extended  In  this  section  to  the  cases 
of  seml-lnflnlte  and  finite  layers  of  stratified  fluid.  The  extensions 
are  based  on  the  normal  mode  expansions  for  these  cases  that  were  derived 
in  [12,  §9].  Only  the  principal  concepts  and  results  are  formulated 
here  since  the  proofs  are  entirely  analogous  to  those  of  the  preceding 
sections. 

Semi-Infinite  Layers.  As  in  [12,  §9]  the  fluid  Is  assumed  to 
occupy  the  domain  and  to  satisfy  the  Dlrlchlet  or  Neumann  boundary 
condition.  Here  the  functions  p(y}  and  c(y)  are  assumed  to  be 
Lebesgue  measurable  and  satisfy 

(6.1)  0  <  <  p(y)  <  Pjj  <  0  <  Cjjj  <  c(y)  <  Cjj  <  " 

and 

(6.2)  |p(y)  -  p(«)|  <  C  y"^,lc(y)  -  c(®)  |  <  C  y“® 

for  all  y  >  0  where  Pj^j,  Pj^,  p(®),  0^^^,  c^j,  c(®),  C  and  a  are  constants  and 

(6.3)  a  >  3/2, 

As  In  [12,  §9]  the  acoustic  propagators  for  p(y)  and  c(y)  corresponding 
to  the  Dlrlchlet  and  Neumann  conditions  will  be  denoted  by  a"  and  A*, 
respectively.  They  are  selfadjolnt  non-negative  linear  operators  In 
JC^  -  L2(R^,c“*(y)  p"*(y)dxdy). 

The  normal  mode  functions  (x,y,p,X)  for  A^,  as  defined  In 
[12,  §9],  are  parameterized  by  (p,X)  e  ■  {(PiX)  [  X  >  c*(*)jp|*}. 
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Their  asymptotic  form  for  y  -►  ®  Is  given  by 


(6.4) 


x,y,p,q) 


c(Ip1»^)  Li(p*x-yq)  +  Rj  gi(p*x+yq) 


2iT 


where  q  -  q(lp|,X)  -  (X  c"*(«)  -  lp|*)*^^  ,  c(1p1,X)  -  (p(»)/4Trq(  |p|  ,X))  , 

■  R^((p1,X)  and  |R^(|p|»X)|  -  1.  As  In  the  preceding  sections.  It 
will  be  convenient  to  Introduce  new  parameters:  (p,X)  (p,q) 

"  (p*q(|p|>^))  €  R^  and  normal  mode  functions 

(6*5)  <l'+(x,y,p,q)  -  (2q)^^*c(®)  (s.y.p.X) 

where 


(6.6)  X  -  X(p,q)  -  c*(®)(|p|^  +  q*). 

The  asymptotic  form  of  Is 


(6.7)  <ji^(x,y,p,q) 


c(”)P  Ll(p*x-qy)  ^  gl(p*x+qy)l 

(27r)  Va  \  J 


y  .*■  oo. 


The  second  family  defined  by 


(6.8)  ♦l(x,y,p,q)  -  <(i^(x,y,-p,q) 


Is  also  needed.  It  satisfies 


(6.9)  ♦_(x,y,p,q)  ~ 


c(«)p^/*  (®) 

(2Tr)V? 


I^Kp'jrtqy)  ^  gl(p«x-qy)J^  ^ 


The  expansion  theorem  of  [12,  §9]  implies  that  the  limits 

M  .  - 

(6.10)  f^(p,q)  •  L2(Rh-llm  ♦+(x,y,p,q)  f(x,y)  c“*(y)  p“‘(y)dxdy 


»f*a> 


exist.  Moreover,  if 


44 


(6.16)  v-(t,x,y)  -  <l>l(x,y,p,q)  exp  (-it  w(p,q))  h^(p,q)dpdq 

Jp » 


where  u>(p,q)  -  c(“)/|p|^  +  q^.  Moreover,  as  in  §3  one  can  write 


(6.17)  <l>^(x,y,p,q) 


c(«»)pl/2  (00) 


(2ir) 


3/2 


I^Kp-wqy)  iJ(y,p,,,+e«i'-*-'»')R 


where 


(6.18) 


lim  I^(y,p,q)  -  1, 

jr^oo 


lim  R^(y,p,q)  -  R^(p,X). 

Then,  proceeding  as  in  §3,  one  can  prove  the  following  analogue  of 
Theorem  3.1. 

Theorem  6.1.  For  every  h  6  let  v^(t,*)  be  defined  by 
(6.19)  Vj(t,x,y)  -  exp  (-itc(<»)  aJ^*)  hp(x,y),  (x,y)  e  R^ 

where  ho  €  L2 (R^)  is  the  function  whose  Fourier  transform  is 


(6.20) 


Then 


ho(P»<l)  ■  1 


c(«)  y(p,q),  (p,q)  e  r|, 

0  ,  (p,q)  €  R^. 


(6.21) 


11m  |Vj(t,»)  -  v|(t,0|„  •  0. 
t"***  ”+ 


^(y.p.q) 
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Transient  Guided  Waves.  For  both  semi- Infinite  and  finite 
layers  the  form  of  the  guided  components  is  precisely  the  same 

as  for  the  case  of  an  Infinite  layer.  Thus  the  analysis  of  §4  applies 
unchanged  to  these  cases. 
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